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I investigate the properties of spontaneous parametric down-conversion (SPDC) involving co- linear 
Gaussian spatial modes for the pump and the photon collection optics. Approximate analytical and 
numerical results are obtained for the peak spectral density, photon bandwidth, pair collection 
probability, heralding ratio, and spectral purity, as a function of crystal length and beam focusing 
parameters. I address the optimization of these properties individually as well as jointly, and find 
focusing conditions that simultaneously bring the pair collection probability, heralding ratio, and 
spectral purity to near-optimal values. These properties are also found to be nearly scale invariant, 
that is, ultimately independent of crystal length. The results obtained here are expected to be useful 
for designing SPDC sources with high performance in multiple categories for the next generation of 
, SPDC applications. 
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I. INTRODUCTION 



Spontaneous parametric down-conversion (SPDC) is a process widely used to produce both entangled photon 
pairs and heralded single photons. SPDC photon sources play a key role in quantum cryptography, optical quantum 
computing, quantum metrology, and in fundamental studies of quantum mechanics. In many applications, a nonlinear 
• optical crystal is pumped by a focused laser beam and the emitted SPDC light is collected into a pair of optical 
^ \ waveguides, such as an integrated optical circuit [I| or a fiber network [21. Efforts over the past decade have led to 
O^. some very bright SPDC sources 0, |j| . As brightness has improved, interest has progressed from simply making SPDC 
sources brighter to controlling other properties of the emitted photons, such as the spectral entanglement in each 
, photon pair [^-[1! ■ While several papers in recent years ^9nl6i] have addressed the question of how to focus the pump 
^ ' and/or collection optics optimally, some important questions remain. Since different studies have invoked slightly 
, different assumptions and optimized slightly different measures, it is not clear how pump and/or collection focusing 
simultaneously affects all the quantities that are generally of interest, and what trade-offs (if any) exist between these 
quantities. Furthermore, the scaling laws for the optimized quantities and optimal parameter values are not readily 
apparent (e.g. whether it is possible to halve the crystal length and, through a change of beam parameters, obtain 
' similar performance). 

To address such questions, I present here a new study of SPDC for the case in which the pump and collecting 
optics define co-linear Gaussian spatial modes. While the theory is general, the envisioned context is SPDC in a 
periodically poled nonlinear crystal with emission in the visible or telecommunication spectral range (a common and 
useful configuration). In this study I consider five properties of the collected biphoton state that are commonly 
of interest: the joint spectral density; the photon bandwidths; the pair collection probability; the heralding ratio 
(pair/single photon collection ratio); and the spectral entanglement. These properties are calculated — in some cases 
d ' analytically, in some cases numerically — as functions of experimental parameters, yielding predictions for the absolute 
values of the properties as well as for the parameter values that optimize each property. Additionally, this study 
reveals several scaling laws and shows that some properties can be jointly optimized while others require a trade-off. 
Many of the results presented here are new, with a few appearing to differ from predictions by others. 

The paper is organized as follows: Section |lT] establishes the foundational equations describing Gaussian optical 
modes and the quantum physics of SPDC. Sections IIIIIIVIII derive expressions for the five properties mentioned above 
and discuss their dependence on experimental parameters. Section IVlIII discusses the results of this study in light of 
prior works, and Section [IXI summarizes the main conclusions. 
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II. SPDC WITH GAUSSIAN SPATIAL MODES 



SPDC is the lowest-order effect of parametric interaction between a strong pump (p) field and two other fields, 
designated as signal (s) and idler (i), initially in the vacuum state. The quantum Hamiltonian governing this interaction 
is 

H = y"eo (x^^Hr) : E+(r,i)E-(r,i)E-(r,i) +H.C.) d^j. (1) 

where r = {x,y,z) is the spatial coordinate, t is time, ^^^•'(r) is the nonlinear susceptibility tensor, and E+(r,^) ~ 
E^{r,t)^ is the positive-frequency part of the electric field quantum operator. In canonical treatments, the field is 
expanded in terms of an orthonormal set of modes. Since the interest here is in SPDC involving Gaussian spatial 
modes, I choose a modal expansion containing a Gaussian mode for each frequency w. It will be sufficient to consider 
just these Gaussian modes until Section IVlII In this case E+(r,i) may be written as 

/•oo 

E^(r,i)= / duj\ Etj(r)e^*'^*atj(t) -I- irrelevant non-Gaussian modes (2) 

Jo V 47rc£o 

where E^(r) is the electric field function for the (Gaussian) mode of frequency w and ajj(i) is the operator that creates 
a photon in that mode. 

The state resulting from the interaction ([T]) is obtained by applying the operator exp — ji dt H (t) to the 

initial state. The part of the state corresponding to SPDC, that is, the creation of a single pair of photons, is just the 
first-order term: 



|*SPDC>=~^/ dii/(t) I initial). (3) 

Putting (P) and © into ^ gives 

l^'sPDc) = -* / duJsdujiij{uJs,uji)a'l^dl.\va,c) (4) 
Jo 

where |vac) is the vacuum state of the signal and idler and 



oo 



■0(t^s,Wi) = ^ ^^^-^^s(wp)C'(a;s,a;i) (5) 

is the SPDC amplitude. Here Eq is the vacuum permittivity, Xj — lixcjiOj is the free space wavelength of field j (j' — 
p,s,i), and Wp = + cjj. The pump field has been assumed to be in a coherent state with independent spectral and 
spatial dependence. Accordingly, the operator a^^ has been replaced by the complex amplitude a(wp) = s{ujp)^y N^, 

where s(wp) is the pump spectral amplitude (with J |s(a;)|^ dui = 1) and Np is the mean number of pump photons 
(the pump energy divided by huip). The quantity 

0{us,u,)^ J x<'>(r):E^^(r)E:^(r)E:^(r)d3r (6) 



medium 



is the spatial overlap of the pump, signal, and idler modes in the medium, which generalizes the sine phase matching 
function encountered in plane- wave treatments of SPDC [l^ . The medium is taken to be a bulk material of length 
L centered at the origin, with cross section large enough to contain the mode functions. The material may be 
ferroelectrically poled so that X*-^^(r) alternates sign with spatial period A. 

The properties of the SPDC state can be calculated once particular forms are chosen for the pump spectrum s{iOp) 
and the mode functions E;^^. (r) . Consideration will be restricted to modes of the form 



h IKZ 



(7) 



which describes a linearly polarized, paraxial Gaussian beam with waist at the origin, propagating along the z axis in 
an optically uniform medium. Here w is the waist size, e is the polarization unit vector, k = nuj/c is the wavenumber. 
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n is the refractive index, and q = w'^ + 2iz/k. This choice is in part motivated by the fact that some of the best SPDC 
sources in existence involve Gaussian beams co-propagating along one of the principle refractive axes of a transparent 
crystalline material, such as potassium titanyl phosphate or lithium niobate. (The anisotropy of the refractive index 
may be safely ignored in such cases.) Additionally, symmetry indicates that the spatial overlap is largest when the 
modes co-propagate and have their waists co- located at the center of the crystal. This leaves the waist sizes Wp,Ws, Wi 
as the free parameters for mode optimization. 



Using modes of the form ([7]), the mode overlap ^ is 



(x^ + y^) ( — + — + — I + i{Ak + mK)z 



dxdydz. (8) 



Here Afc = kp — kg — ki is the wavenumber mismatch, K — 2tt/A is the poling spatial frequency, m is the order 
of quasi-phase matching, x'cs = x'^' • GpegGi is the effective nonlinear coefficient, and e is an efficiency factor that 
includes Fresnel loss and the Fourier coefficient of the mth harmonic of the poling spatial function. Integrating over 
X and y yields 

Before proceeding, it will be convenient to rewrite (|9]) in terms of dimensionless quantities. The primary independent 
variables are the phase mismatch 

^={Ak + mK)L (10) 

and the focal parameters 

where £_j 1 (^ 1) means that field j (j =p,s,i ) is focused strongly (weakly) relative to the length of the crystal. 
I also define the auxiliary quantities 

A,.l + ^|i±|L|L, (12) 

"-p Sp "-p Sp 



Afc Ah 



C'^—i^^. (14) 



and the aggregate focal parameter 



^+ ?p 



The quantities j4+, and ^ are independent and uniquely determine ^s, ^i, and ^p. C is determined by A-^- and 
B+. Note that A±, B±, and C do not depend on the absolute values of the focal parameters, but on the focus of the 
signal and idler relative to the pump. In terms of these dimensionless quantities, eq. ([S]) may be written as 



^(.s,.0 = /^--^^--sr.iVpLx| ^^^P^dl. (16) 



eoup AsAi y/A+B+ 7-i 1 - - CS,"^! 

Eq. will be the starting point for analysis in sections IIIHIVIl 

Four approximations will be invoked throughout this work in order to simplify analysis and yield more useful 
results. These approximations are reasonable for typical bulk SPDC sources, which may be defined to have the 
following characteristics: the length of the medium is > 1mm and its refractive index is > 1.5; the parameteric 
interaction is quasi-phase matched with a first-order grating of period A > 5 /im; and emission is in the visible or 
telecommunication spectral range, with As ^ 1.6 /xm and Ap < 0.8 /um. When assessing the accuracy of approximate 
formulas, these values will be considered to represent the worst typical case. I will also consider a "reference source" 
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consisting of degenerate type II SPDC in 10 mm periodically poled potassiam titanyl phosphyate (PPKTP) with a 
750 nm pump. This is similar to several sources that have been demonstrated to have good performance [sl, R fisl . 
One approximation that will be used liberally is 

,17, 

This approximation is motivated by the fact that efficient SPDC occurs when Afc « —mK, where mK/kj is typically 
much smaller than 1. 

Another approximation that will prove convenient is 

CwG. (18) 

The actual value of C depends in a complicated way upon the experimental parameters, but can be shown to obey 
the bound 



\C\< 



|Afc| fcp 
Aksh 



AkshL iksh 



(19) 



when (|17l) is valid. One then has \C\ ^0.1 near phase matching conditions, which turns out to be small enough to 
make (|18p a fair approximation over the range of focusing that is desirable with respect to the five properties of the 
state considered here. 

The third approximation is that the frequency dependence of (jl6p is determined essentially by the pump spectrum 
s{ujp) and the frequency dependence of the wave mismatch Ak. That is, I take 

(2) 



^±,S±,Om/™T^~ constant (20) 

over the range of (cjs,i^i) for which the amplitude is appreciable. For (|20p to hold, the bandwidths of the photons 
must be much smaller than an optical frequency. Type II and non-degenerate type I SPDC sources typically satisfy 
this condition, but frequency-degenerate type I sources and sources made of very short crystals, which tend to have 
large bandwidths, may not. 

Finally, after section lllll it will be assumed the wave mismatch Afc has a predominantly linear frequency dependence: 

6kj « -^Sio^ (21) 

where n'j = cdkj/duj is the group index of mode j and Sluj {5k j) denotes a shift from the nominal frequency 
(wavenumber) of mode j. Again, this approximation is generally valid for type II and non-degenerate type I SPDC 
sources, but not for frequency-degenerate type I sources. 

The impact of these approximations, particularly (|18p . will be addressed more fully in the context of each major 
result. 



III. JOINT SPECTRAL DENSITY 

The joint spectral density |7/;(a;s, Wi)|^ is the expected number of photons pairs, per signal bandwidth per idler 
bandwidth, emitted into the Gaussian collection modes. If the collected photons pass through spectral filters of 
narrow bandwidth, the effective brightness of the source is determined by the joint spectral density at the filter 
frequencies. Let us now consider the problem of maximizing this quantity, that is, finding the values of ^p, ^s, 6j ^ that 
maximize \ij){uJs,to\)f' given the frequencies Wp, Wg, i^i and crystal length L. From p6p . it is apparent that maximization 
determines certain values (^max^ ^max^ ^max^ (j^max^^ The beam parameters Cp,^s,6 are then complicated functions 
of (/cp, fcs, fci, ^+'*'^, 5+'^'', ^'™^'') via relations (fT0|) - (fT5)) . However, substantially more informative results can be 
obtained if the P term in is neglected. With C w (approximation (IT51) '). the function to be maximized is 

1 r VeexpW2]^^_ ^22) 



The two factors are independent and can be maximized separately. A^B^ can be written as 

= (^1 - (1 + A,r, + X,n) (^1 + ^ + |i 
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where 



kj / 1 + Ak/kj 
k^V l~Ak/kp' 



/ I- Afc/fcp 

"■^^ - V 1 + ^k/kj ■ 



— 1 /2 

The maximum of occurs at = n = 1- Invoking approximation f|17p gives 



6 « 6 « ?p « ? 



and 



max ■ 



The second factor in ([^ . 



J-1 1 - 



(23) 
(24) 

(25) 
(26) 

(27) 



requires numerical evaluation and is plotted in Fig. [T] (Note that for ^ ^ 1, is proportional to the usual phase 
matching function sinc($/2).) The maximum value of F is 2.06 • • • at ^ = 2.84 • • • , $ = —(1.04 • • • )7r. This result, 
together with (|26p . yields the maximum spectral amplitude 



max |V'(ws, Wi)| « 1.03< 



eoUp AsAi 



V^s(wp) 



which is obtained under the conditions 



6 « Ci « Cp 



2.84, 



$ w -1.047r. 



(28) 



(29) 
(30) 



It should be noted that these are precisely the conditions which have long been known to maximize the efficiency of 
sum frequency generation and parametric amplification with Gaussian beams p^ . The spectral density |^/'(a;s, "^i)!^ 
is proportional to the number of pump photons A^p and to the crystal length. (In contrast, for the case of collimated 
(plane wave) interaction, the spectral density would grow quadratically with the crystal length.) It may also be noted 
that, for all other factors being equal, short-wavelength sources are brighter than long-wavelength sources. 

Since it may not always be possible to implement or verify condition (|29p accurately, it is worth examining how 
sensitive the peak spectral density is to the focal parameters. Specifically, I compare max$ 1-01 , which is a function 
of (Cp,^s,6)j to the global maximum value max^_$ |-0| . I constrain the parameter space by the condition = ^i, 
which is not only optimal in regard to the spectral density, but as will be shown in later sections, is also optimal or 
near-optimal for other quantities of interest. Fig.[2]shows the relative brightness max$ /max^_$ 1-01^ as a function 
of and — It can be seen that the regime of good focusing is rather broad: the waist size (and correspondingly, 
the angular divergence) must be made smaller or larger by roughly a factor of 5 to reduce the peak spectral density 
to half its maximum value. Over this range the optimal phase $ (not shown) varies from —tt/2 to — 37r/2, meaning 
that the frequency of the peak shifts slightly with focusing. If the signal and idler focus are fixed, the optimal pump 
focus is given by ^p = Cs(= 6)- However, if the pump focus is fixed, the optimal signal and idler focus are given by 

= 6 « v/2:8ie^. 

These results are robust to approximation: for the worst case defined in section |lTl the "optimal" values given in 
([29]) and (f30|) differ from the true optimal values by 5% and 10% respectively, but, due to the broadness of as a 
function of these parameters, yield an amplitude that is within 1% of the actual maximum value. 



IV. PHOTON BANDWIDTH WITH MONOCHROMATIC PUMPING 



SPDC photons are sometimes employed in systems with limited optical bandwidth. In designing sources for such 
systems it is helpful to know how the photon bandwidths depend on the design parameters. One case of particular 
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FIG. 1: (color online) The spatial overlap factor _F'(^, $) defined in eq. (I27p . The cross marks the location of the peak. 
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FIG. 2: (color online) Peak joint spectral density, normalized to the global maximum, as a function of the pump and photon 
mode focus. 

interest is when the pump bandwidth is chosen to maximize the spectral purity of the photons. That case will be 
addressed in section IVII Another particularly interesting case is that of monochromatic pumping, which will be 
addressed now. 

With a monochromatic pump, and under approximation ()20p . the spectral dependence of ()16|) arises primarily from 
the dispersion of the phase mismatch $ = (Afc + mK)L. The phase mismatch may be expanded as Taylor series in 
5u)s and 5uJi, the deviation of the signal and idler frequencies from their nominal values. Assumption (j2ip then gives 



$ w <i>o 



— {5ujs + <^i^i) -5los 

c c 



(31) 



where $o is the phase mismatch at the nominal frequencies. Monochromatic pumping constrains the frequencies to 
Suis = — yielding 



$ w $0 + —Sujs 



(32) 



The fact that <i> is (approximately) a linear function of ujs means that the photon bandwidth Awg {— Awj) can be 
expressed in terms of a dimensionless "phase mismatch bandwidth" A$, which is the width of when expressed 
as function of $. 
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FIG. 3: (color online) The normalized photon bandwidth as a function of focusing. The dashed line is the heuristic formula 
(133}. 



Under approximation C « 0, is proportional to |F((^, <i>)p. The full-width at half maximum (FWHM) of 



\F{£^, $)|^ (see Fig.[T|) was calculated numerically and is plotted as the solid line in Fig.[3l I find that the behavior of 
A<i> in this case is captured well by the heuristic formula 

A$ - 27rmax(l,^/10), (33) 

shown as the dashed line. Combining ([32]) and (l33l) yields the photon bandwidths 



where b = L/£_ is the aggregate confocal length of the modes (reducing to fcpWp when = = ^p). When focusing 
is weak to moderate (^ < 10), the photon bandwidth is determined by the crystal length (the well-known 1/L 
dependence). But when the focusing is strong (^ > 10), the bandwidth is larger [20^ and determined instead by the 
confocal length, going as 1/5. 

The reason that tightly focusing the pump increases the bandwidth can be understood as follows: At negative 
values of the phase mismatch (i.e . away from the nominal wavelengths), the photon spatial distribution in the far 
field takes the form of a ring [21]. With a collimated pump, a certain change in wavelength makes the ring radius 
larger than that of the collection modes, and photons of those wavelengths are not collected. But when the pump 



is tightly focused, the rings are spatially broadened and the photons partially overlap the collection mode [22|. The 
larger the spatial broadening, the larger the wavelength change must be for the photon distribution to lie outside the 
collection mode. 

Due to the approximation C = 0, the bandwidth plotted in Fig. [3] is not exact. For the reference source, the error 
in Aws is < 3% for ^ < 10 and < 10% over the range 10 < < 100. For the worst typical case, the error is < 3% 
for ^ < 1, but increases to 10% at ^ = 10 and varies between 30% and 50% over the range 10 < ^p < 100. Thus 
bandwidth predictions may not be highly accurate in the regime of focus-induced spectral broadening. Nevertheless, 
even in the worst case the approximate SPDC state (calculated with C = 0) has a large overlap with the exact 
state l^/i) (calculated using the fidehty (V^jV') V(^I^)(V'IV') is 0.999 at ^ = 1, 0.97 at ^ = 10, and 0.91 at ^ = 100. 

Also, the heuristic formula A$ ^ 27rmax(l, ^/^o) remains as accurate as shown in Fig. Ofor a suitably chosen value 
of Co. 

V. PAIR COLLECTION PROBABILITY 

When SPDC photons are collected without spectral filtering, maximizing the pair collection probability 

^si = / \ip{uJs,i^i)\'^ dujsduji. (35) 
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FIG. 4: (color online) Dependence of the pair collection probability on the aggregate focus. The dashed line is eq. (|40p 
normalized by its asymptotic maximum. The solid lines are calculated using "worst case" values for the parameter C (see 
discussion surrounding eq. H19|) ) and normalized by the same factor. 



is usually more important than maximizing the peak spectral density. (In most experiments, Psi is proportional to 
the coincident photodection rate, or the probability of detecting a pair of photons following any given pump pulse.) 
Under approximation (j20p . the dominant spectral dependence arises from s{LOp) and Afc. The phase mismatch (|3ip 
may be written as 



2c 



where w_ 



Wj. Then diUscLji 

r 



{duj^dujp)/2 = d^dujp c/{\n'^ 



\<-n[\L 
STr^hcensUi 



/ (2) \ ^ 



n 



s(tjp) 



2c 

"il-^): giving 



dl 



dujpd^. 



Under the approximation C ~ 0, the only place $ appears is in the exponential function. 
/ exp - l')/2] d$ = 47r5(/ - /') and / |s(a;p)| dujp ^ 1, one obtains 



Psi 



32TT^hcensni 



which may be directly integrated to give 



AgAi 



dl 



+ J-i 



Psi 



AgAi 



arctan (^) 
A+B+ 



(36) 

(37) 

(38) 
Using 

(39) 
(40) 



Formula ([Ml. plotted as the dashed line in Fig. 21 suggests that the pair probability can only be optimized in an 
asymptotic sense — that there are no finite values of , , 6 that maximize Pgi . In reality, the asymptotic behavior 
holds only as long as £,kp / (ALksh) < 0.1. With very tight focusing or very short crystals, the approximation C — 
breaks down, causing Psi to peak near its asymptotic value (solid lines in Fig. U). In the worst typical case the pair 

probability error 1 - {^j^lip) / is < 4% for ^ < I, < 13% for ^ < 10, and < 20% for ^ < 100. For the reference 

source described above, the pair probability error is < 1.5% for ^ < 100. In any case, Pgi evidently has an upper 
bound 



Psi < 



Sir^hcensrii 



' (2) ' 
Xcff 



EQUp \n'^ -n[\\ AgAi 



Nr. 



(41) 
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FIG. 5: (color online) Dependence of the pair collection probability on the focus of the pump and photon modes. 



that cannot be exceeded, no matter how long the crystal. (This somewhat surprising result will be discussed in section 
IVIIin For type II SPDC sources made of PPKTP or periodically poled lithium niobate (PPLN), eq. (gT]) predicts 
that brightnesses exceeding 10~^ collected pairs/pump photon should be achievable — roughly an order of magnitude 
brighter than the brightest existing sources. 

The dependence of Pgi on and = 6 is shown in Fig. [5] If the signal and idler focus are fixed, the optimal 
pump focus is given by = ^s(= 6)- If instead is fixed, the optimal value of is nearly equal to for ^p > 10 
and slightly larger than ^p for ^p < 10. 



VI. SPECTRAL PURITY 



A number of applications of SPDC photons, including the generation of multiphoton entangled states for quantum 
computing |23l426l |. involve interference between photons from separate (and nominally identical) SPDC sources. The 
success of such applications depends not only on the efficie ncy of SPDC photon production, but also on the degree 
of mutual coherence of photons from independent sources [23, . This coherence, which is directly related to the 
interference visibility, is given by the single-photon purity 

P = (42) 



where 



-0(^8, Wi) = y2y^Uj{uJs)Vj{uJi) (43) 



j 

is the Schmidt decomposition [2^ of the collected biphoton state. The single-photon purity is inversely related to the 
degree of entanglement between the signal and idler frequencies, with p = I corresponding to no spectral entanglement. 
(Note that any spatial entanglement present in the emitted state is discarded upon post-selecting photon pairs that 
couple into the Gaussian collection modes.) 

It is natural to ask whether the parameters that yield high source brightness also yield high spectral purity. Under 
approximations ([20|) and ([21]), the frequency-dependent part of ip{u}s,u}i) is s(wp)i^(^, $). Since the purity depends 
only on the shape of ip{ujs,uji), and not on its location or extent within the (a;s,a;i) plane, Wg and may be replaced 
with the dimensionless variables ujs/fl and uji/fl, where 57 = c/{n'^ — n[)L. Given a particular functional form for the 
pump spectrum, the purity may then be computed as a function of (^, 9^ Aup/fl) where Alo is the pump bandwidth 
and 

tan9=— ^ (44) 
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FIG. 6: (color online) Dependence of the spectral purity of the collected photons on the mode focus. The numbers accompanying 
each line are the group velocity angle 9 defined in eq. (|44|l . Solutions for 6 > 45° mirror those for 9 < 45°. 

is the slope, in the (ws,a;i) plane, of the line characterized by $ = 0. Note that since ^ is the only relevant focusing 
parameter, we are free to set ~ Ci ~ to maximize the brightness at a given {^,9, Acop/fl). Also, note that the 
brightness is independent of pump bandwidth when approximation (j2ip is valid. 

The spectral purity obtainable with a Gaussian pump was determined by computing s{ujp)F{^, $) over a grid of 
sufhcient extent and resolution in the {iOs/^,uji/fl) plane for various parameter values {^,0, Aujp/n). For each set of 
parameter values, singular value decomposition was applied to the matrix of computed values to determine {ffj} and 
p. Numerical optimization was then used to determine the value of Aujp/D, that maximizes p at each {^,,0). In Fig. [6] 
these optimized values of p are plotted as a function of ^ (= Cp) for several different values of 9. For 9 = 45°, the peak 
purity is 0.94 and occurs at ^ = 2.2. As 9 decreases to 0°, p asymptotically increases to 1 for all values of ^, although 
the optimal pump bandwidth becomes infinite in this limit. Below 0°, p drops rapidly. Solutions for 9 > 45° mirror 
those for 9 < 45°, with the signal and idler exchanging roles. 

These results can be understood as follows. The purity is directly related to the factorability of tp{ujs,(j-'i) oc 
s(cjp)i^(^, $) into separate functions of and Wi. When the pump function s(wp) is Gaussian, complete factorability 
can be achieved if and only if the phase matching function F{^, $) is also Gaussian with orientation < < 7r/2 p. 
But F{^, $) is not Gaussian: for ^ <C 1 it is a sine function, which has oscillating side lobes, while for f ^ 1 it has 
large skew. In both these regimes, the purity is reduced. F{^, $) is most nearly Gaussian for ^ ~ 2, which becomes 
the regime of greatest purity. 

Comparing with Figs. [5] and SI one sees that in the regime of peak purity ^ 2), the spectral intensity is close 
to its peak value and the total collection probability is more than 70% of its asymptotic maximum (PT|) . Thus, high 
spectral purity (> 94%) can be obtained with focusing conditions that also yield relatively high brightness, provided 
< 9 < 90° (which amounts to the requirement that rip lie between n'^ and n[). 

An explicit evaluation of the accuracy of these results in view of approximation (jlSp was not performed. However, 
it may be noted that the results of this section are based on a decomposition of an approximate state which, as 
discussed at the end of section IIVI has a high degree of overlap with the actual state. 

VII. SINGLE-PHOTON COLLECTION AND HERALDING RATIO 

In some applications of SPDC, the detection of a photon in the signal mode is used to indicate the (probable) 
presence of a photon in the idler mode. (Although the photons are always emitted in pairs, conditions generally 
allow one photon to be emitted into a spatial mode defined by the collection optics while its partner is emitted into a 
non-collected spatial mode.) In such applications it is desirable to have a high heralding ratio % = Psi/Ps, where Pg 
is the probability that a photon is emitted into the signal collection mode, regardless of the mode its partner is in. 
Since the probability of collecting both photons cannot exceed the probability of collecting one of the photons, the 
heralding ratio can be at most unity. In applications that treat signal and idler photons symmetrically, the quantity 
rjsi = Psi/y/PsPi is often the metric of choice. In this work rjs will be called the signal heralding ratio and rjsi will be 
called the symmetric heralding ratio. 
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To obtain the signal collection probability Pg, the joint collection probability must be summed over a complete set 
of spatial modes for the idler photon. The Laguerre-Gauss modes 



(r;w) = 



V2 



i+i 



2 ^2 



2wfp 



exp 



h ikiZ + il(f> 

1i 



(45) 



form such a set, where is the {n,l) associated Laguerre polynomial, p — ^ +2/^, and tan(/) = xjy. The 
fundamental mode with n = / = is just the Gaussian mode introduced in section |lTl Since the pump and signal 
mode are azimuthally symmetric, the spatial overlap vanishes unless the idler mode is also azimuthally symmetric 
{I = 0). Thus the spatial overlap involving the nth idler mode is 
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With the Laguerre polynomial expansion formula L„(x) — X]j=o(~-'^)"' 
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tX^ and a bit of work, one can obtain 
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Applying definitions (|10l) - p^ and making the approximation C « yields 
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In analogy with formulas psp and ([5|), the signal photon probability may be written as Pg = 

/ IV'ti(ws, Wi)|^ dojsduji where ■;/'„ = y/2TT^hNp/eoXpXsXis(ujp)On{(j-'s,i^i)- Following the approach taken in section 

IVl we have 
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Integration yields the signal collection probability 
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(53) 



The corresponding formula for the idler probability Pi can be obtained by interchanging the labels s and i everywhere. 
Eq. (|5^ is very similar to eq. (HUl) and holds under essentially the same conditions. Like Psi, the signal probability 
Ps is (to first approximation) an asymptotically increasing function of focal parameters, with an upper bound 



32Tr'^hcensni I Xpg 
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Ps < 



Np. 



(54) 
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FIG. 7: (color online) Dependence of the signal collection probability on the focus of the pump and signal modes for near- 
degenerate SPDC (fcs ^ki^ kp/2). 



Also like Pgi, Ps is locally maximized by taking w (see Fig. [5]); however Pg varies more slowly with and 
than Psi (Fig. ^ . The fact that Ps is broader than Pgi means that collection of the signal in a non-optimal Gaussian 
mode projects the idler onto a mode that does not couple well to a Gaussian mode of any size, i.e. a mode that is not 
Gaussian. Note also that (l53l) is independent of as it should be (the parameters of the idler collection mode should 
be irrelevant when one does not care whether the idler is collected). 

If one takes = 6 = Cpi which locally maximizes both P^i and Pg, the signal heralding ratio reduces to the simple 
expression 

An analogous expression exists for rji. These expressions give rjs — rji = 0.75 for near-degenerate SPDC (for which 
fcg sa fci kp/2) and values less than 0.75 for the non-degenerate case. However, higher heralding ratios are possible 
with different focusing conditions. The optimal source configuration with regard to both heralding and brightness 
is not a single set of parameter values, but a curve in parameter space having the property that ijs (or rjsi, if that 
is of interest) cannot be increased without decreasing Pgi, and vice versa. Numerical methods were used to find the 
values (^s, ^ii ^p) that maximize either rjs or ?7si at 50 values of Pgi covering the range from to maxPgi. The resulting 
points are plotted versus ^p in Fig. |8] for the case fcg ~ M ~ fcp/2. Fig. [8] shows that it is possible to achieve very high 
heralding ratios, but only with a substantial reduction in brightness: a factor of at least 4 to achieve % > 0.95 and a 
factor of at least 10 to achieve ?7si > 0.95. That is, there is a trade-off between brightness and heralding ratio, which is 
somehwat worse for symmetric heralding than asymmetric heralding. Both r]s and 77si approach unity in the limit that 
the pump is collimated (i^p — > 0). In this limit the best trade-off between Pgi and rjsi is achieved with = 6 ^ Cpj 
while the best trade-off between Pgi and 77s is achieved with ^ 3^s ~ 3fp. The trade-offs between Pgi and 77s or r]si 
are found to be slightly worse in the non-degenerate case than the near-degenerate case. 



VIII. DISCUSSION 



Perhaps the most significant finding of this study is that there exist focusing conditions (Cp,Cs,6 ~ 2.5) that 
simultaneously bring the brightness, heralding ratio, and spectral purity to substantial fractions of their maximal 
values. Whenever n'p lies between and n[, a spectral purity of at least 0.94 can be obtained with at least 74% of 
the maximum achievable brightness, while the heralding ratio is at least 0.75 (for near-degenerate SPDC). Maximum 
brightness can be achieved with tighter focusing at the cost of reducing the spectral purity (to 0.7 at worst) The 
worst trade-off is between brightness and heralding ratio: focusing conditions that yield a symmetric (asymmetric) 
heralding ratio of 0.95 or better reduce the brightness to 10% (25%) of its maximal value; however the spectral purity 
remains 0.84 or better. 

Another notable finding of this study, which to my knowledge has not been reported previously, is that many 
important properties of the collected SPDC state are essentially scale invariant. To the extent that the parameter C 
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pump focus pump focus 

FIG. 8: (color online) (a, 6) Simultaneous optimization of the total collection probability Psi and the symmetric heralding ratio 
r^ai. (c,rf) Simultaneous optimization of Psi and the signal heralding ratio rjs- Panels (b) and (rf) plot the focal parameters that 
yield the optimal curves shown in (a) and (c), respectively. In (a) and (c) the collection probability goes with the left axis and 
the heralding ratio goes with the right axis. 

(see (|14l) ) is negligibly small, the joint collection probability, heralding ratios, and spectral purity are determined by 
the dimensionless ratios L/kjw"^ and AwpL/c. Changing the crystal length has no effect on these three quantities if 
the pump duration and confocal ranges are also scaled by the same factor. In the end, the crystal length merely sets 
the bandwidth of the system, with longer crystals yielding (or requiring) smaller bandwidths. The idea that a longer 
crystal does not yield a brighter source is perhaps surprising, since longer crystals yield higher SPDC efficiencies when 
the pump is coUimated or weakly focused. Why this is not so for optimally focused sources may be understood by 
noting that, although a shorter crystal has a shorter interaction length (which decreases the spatial mode overlap), it 
also allows the modes to be focused more tightly (which increases the mode overlap). Of course, this argument does 
not hold for arbitrarily short crystals; at some point, the focusing becomes so strong and the bandwidths become so 
large that the paraxial approximation implicit in (jT]) , as well as approximations such as (j20l) , break down. At the other 
extreme, very long crystals may also show worse than expected performance due to the challenge of manufacturing 
very long crystals of high quality. 

Some of the results obtained here appear to agree with prior works, while others appear to differ from prior works. 
The focusing condition = 6 = = 2.84 which maximizes the joint spectral density of SPDC is the same as that 
found by Boyd and Kleinman for maximizing second harmonic generation [l9j . This makes sense, since SHG with 
a monochromatic Gaussian input beam produces a monochromatic Gaussian second harmonic field; optimization of 
SHG then amounts to finding the parameters of the Gaussian modes that maximize the spatial overlap, which is 
precisely what was done in section Hill More recently, Ljunggren and Tengner [13 ] have performed numerical studies 
yielding the optimizing conditions 2^p = 1.7, 2^s — 2^; — 2.3, and the prediction that the optimized joint probability 
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goes as \fL (rather than being independent of L as claimed here). These discrepencies may be due to two key 
differences in approach: In |l4| . optimal focusing conditions are obtained for $ = 0, whereas here and in [l^ the 
condition $ = — 1.047r is found to yield a slightly higher spectral density. Additionally, [l3| employs a plane- wave 
approach in which the diffraction of the pump mode is effectively ignored; this approximation is also made in fis} , 
[ll| and [l5| . In the present study the pump is treated as a diffracting paraxial beam. Experiments to date are 
insufficient to resolve these differences as they cover a limited range of crystal lengths and focal parameters, and are 
complicated by poor repeatability (changing the focus of a beam, or replacing a crystal with one of different length, 
generally necessitates realignment). 

The results of this study are intended to be useful for designing and predicting the approximate performance of 
a promising class of SPDC sources. A number of points should be kept in mind, however. Firstly, all the results 
presented here apply to the collected part of the biphoton state; many features of the emitted SPDC light are lost 
or substantially altered when the state is projected onto the Gaussian collection modes. Secondly, while care was 
taken to confirm the general viability of the approximations employed, it is not difficult to envision sources that would 
violate the assumptions of this analysis and exhibit quantitatively or qualitatively different performance. In particular, 
results obtained here may not be reliable for very short sources [L < mm), sources with very short poling periods 
(A < 5 /xm), and/or sources employing quasi-phasematching of order m > 1. Thirdly, certain physical details deemed 
inconsequential — such as the vector diffraction of the Gaussian beams, the optical anisotropy of nonlinear crystals, 
and the proper normalization of electromagnetic modes quantized in a dielectric — were simply ignored. Fourthly, a 
limitiation of this study is that it does not apply to angle-tuned sources in which beam walk-off plays a role. However, 
since walk-off decreases spatial overlap, it is not clear that such sources could achieve better overall performance than 
the co-linear sources considered here. Finally, of all the results obtained here, only those in section Hill are potentially 
applicable to frequency-degenerate type I sources. The quadratic relationship between frequency and phase mismatch 
in these sources substantially complicates analysis; however, the possibility that these sources might have different 
scaling laws makes them potentially worth a comparable study. 



IX. SUMMARY 



In summary, I have presented here a new theoretical study of SPDC addressing multiple properties of the emitted 
photons that are important in various applications. Analysis was restricted to the promising class of SPDC sources 
involving focused, co-linear Gaussian modes for the pump field and collected photons. Analytical and numerical 
calculations yielded approximate predictions for the peak spectral density, photon bandwidths, absolute pair collection 
probability, heralding ratio, and spectral purity. A scaling law was found which shows most of these properties to 
be ultimately independent of crystal length. It was also found that such sources can simultaneously exhibit high 
brightness (predicted 10~^ pairs/pump photon), high spectral purity (> 0.94), and moderately high heralding ratio 
(> 0.75) when the confocal ranges of the modes are on the order of half the crystal length. Higher heralding ratios 
can be achieved, at the cost of significantly reduced brightness, by focusing the modes less tightly. The results of this 
study are applicable to "typical" SPDC sources, excluding sources of very short length, very short poling period, or 
very large bandwidth. Frequency-degenerate type I sources, which often do not satisfy these criteria, may be amenable 
to a similar kind of study. 
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Directed Research and Development Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for 
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